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Abstract. In this paper we first discuss observational evidence of longitudinal concentrations of 
' magnetic activity in the Sun and rapidly rotating late-type stars with outer convective envelopes. 
Scenarios arising from the idea of rotationally influenced anisotropic convective turbulence be- 
ing the key physical process generating these structures are then presented and discussed - such 
effects include the turbulent dynamo mechanism, negative effective magnetic pressure instability 
■ (NEMPI) and hydrodynamical vortex instability. Finally, we discuss non-axisymmetric stellar 
mean-field dynamo models, the results obtained with them, and compare those with the observa- 
tional information gathered up so far. We also present results from a pure a 2 mean-field dynamo 
model, which show that time-dependent behavior of the dynamo solutions can occur both in 
the form of an azimuthal dynamo wave and/or oscillatory behavior related to the alternating 
0^ energy levels of the active longitudes. 
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The spectroscopic time series obtained for rapidly rotating late-type stars, analyzed 
with Doppler imaging methods, have revealed huge high- latitude temperature anoma- 
lies unevenly distributed over the stellar longitudes. The analysis of photometric time 
t-H ■ series have in many cases confirmed the spectroscopic results, and allowed the follow- 
up of the evolution of the light curve minima for longer time intervals and with better 
time sampling. From these investigations it seems clear that in many stars the magnetic 
activity manifests itself as active longitude(s) that undergo apparently quite irregular 
^vq ■ phase jumps, sometimes called flip-flops, and drifts in the rotational frame of reference. 

In terms of dynamo theory, a change from the solar-type mostly axisymmetric and quite 
regular oscillatory solution to a non-axisymmetric non-stationary mode is expected to 
_ occur as the rotation rate is increasing. 
i-H ■ Some of the phenomena related to these active longitudes can rather straightforwardly 

be understood using mean-field dynamo theory (see e.g. Krause & Radler, 1980). The 
• »"H ■ change from axi- to non-axisymmetric solutions with increasing rotation rate, for in- 
| stance, is a prediction dating back to kinematic linear dynamo solutions, that has been 
later confirmed by more complex non-linear modeling. Time-dependent solutions, with 
the angular frequency of the non-axisymmetric modes being different from the axisym- 
metric mode and stellar rotation constituting an azimuthal dynamo wave, have also been 
recognized as common linearly preferred and non-linearly stable solutions to the dynamo 
equations. Understanding and explaining the flip-flops has remained a prominent chal- 
lenge both observationally and theoretically. Photometry and Doppler imaging naturally 
do not give the information of the polarity of the spots on the active longitudes. There- 
fore it has been impossible to differentiate whether the phase jumps are related only 
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to the alternation of the magnetic energy levels of the spots, or are polarity changes 
also related to the phenomenon. Only very recently such observational data has become 
available; as will be discussed in Sect. 12. 2[ the analysis and interpretation of the data is 
far from straightforward, and clear conclusions are still missing. Theoretically, oscillatory 
solutions are more generally associated with a£l dynamo solutions, where differential ro- 
tation is dominating the field generation over the a effect. In the rapid rotation regime, 
however, differential rotation is expected to be suppressed. It has remained challenging to 
find oscillatory solutions of pure a 2 type, although it has been known since the kinematic 
calculations that such solutions are possible. Recently, more results of oscillatory a 2 dy- 
namo solutions have emerged, and in this paper we also demonstrate one such model in 
Sect. I5XT1 

One of the most puzzling results related to the spectropolarimetric observations, yield- 
ing the surface magnetic field strength and orientation with the surface temperature dis- 
tribution, is that very seldom clear anticorrclation between the magnetic field strength 
and temperature (see e.g. Kochukhov et al. 2012 and references therein), expected from 
the sunspot analogy, is found. From theoretical point of view, it is, indeed, quite poorly 
understood how the dynamo-generated sub-surface fields transform themselves into sun- 
or starspots, and whether the same mechanisms are actually at play in all the objects 
under study. The sunspot sizes, for instance, are on average much smaller than the scale 
of the toroidal magnetic field belt that they emerge according to dynamo theory - this 
is to be contrasted to the starspots that have surfaces areas large enough to be regarded 
as global large-scale phenomena as the dynamo-generated field itself. 

The rather widely accepted paradigm of sun- and starspot formation is the so called 
rising flux tube model (see e.g. Choudhuri & Oilman, 1987), according to which strong 
magnetic flux is generated and stored in the tachoclinic shear layer just beneath the con- 
vectively unstable region. The magnetic field becomes unstable and rises to the surface, 
in the form of thin flux tubes, due to buoyancy force, during which process it becomes 
only weakly affected by the turbulent motions in the convection zone. Closely related to 
this paradigm is the flux transport dynamo scenario (see e.g. Dikpati & Charbonneau, 
1999), furthermore relying on the inductive effect due to sunspot decay near the solar 
surface (Babcock-Leighton effect) accompanied with a conveyor belt provided by merid- 
ional circulation, reviewed by Choudhuri (2013) in these proceedings. Again, convective 
turbulence plays a negligible role in the flux transport model, even the turbulent diffusion 
being reduced significantly from the simple estimates e.g. from mixing length theory. In 
this paper, we concentrate more on the turbulent dynamo picture (see e.g. Brandenburg 
& Subramanian, 2005; Pipin, 2013 in these proceedings), which relies on the idea of mag- 
netic field generation both by rotationally affected anisotropic convective turbulence and 
large-scale flows and non-uniformities in the stellar rotation profile. Sun- and starspot 
formation mechanisms directly related to turbulence effects include the negative effec- 
tive magnetic pressure instability (NEMPI) and hydrodynamic vortex instability, both 
of which are also reviewed in this manuscript. In the context of NEMPI from turbulent 
convection, we also refer to the paper by Kapyla et al. (2013) in these proceedings. 

Next, we give a short review of stellar dynamo models capable of solving for non- 
axisymmetric modes, and the results obtained with them in Sect. 13.41 We note that 
results in this field are much less abundant than in the solar case, where axisymmetric 
(in two dimensions) modeling is adequate. Finally, a summary and main conclusions are 
presented in Sect. |4j 
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2. Observational studies of active longitudes 

The study of longitudinal concentrations of magnetic activity started with the accu- 
mulation of data of the distribution of solar activity tracers on its surface (systematic 
sunspot observations, solar flares, and so on). It was realized that they do not occur 
completely randomly over the longitudes, but 'hot spots' or 'activity nests' surviving 
over several years, maximally over the whole solar cycle, in which the major sunspots 
and flares were accumulating over time, frequently occur. The first tool was to create 
family trees (time-longitude diagrams) of the activity complexes, aiming at following the 
pattern their formed over time on the solar surface - using this technique it soon became 
evident that the hot spots were not always moving with the same rotation rate as the so- 
lar surface, nor the Northern vs Southern activity complexes as a synchronized structure 
(see e.g. Bai, 1987 and references therein). Evidently, something interesting was going 
on, which ignited an intensive study of the phenomenon with sophisticated time series 
analysis tools. 

2.1. Is the solar magnetic field axisymmetric? 

Several parametric and non-parametric time series analysis methods have been used to 
analyze the distribution of sunspots and solar flares over the past few decades. As the 
family tree analysis had already revealed that the rotation period of the hot spots was 
unlikely to match with the surface rotation rate, the first hypothesis was to assume that it 
docs not, and try to find the most suitable rotation period that produces the most uneven 
distribution (i.e. the strongest clustering) of the activity traces over the solar longitudes. 
Again, South and North seemed to be decoupled, the other natural hypothesis being to 
keep them separate in the analysis. Very often the term 'active longitude' is used, in a 
strict sense meaning, that at any latitude, being it Northern or Southern, the activity 
should manifest itself at the same longitude. This definition does not make much sense in 
any observational case - the solar hemispheres behave in a decoupled manner, and usually 
only a fraction of the whole latitudinal extent is seen of other stars. In the dynamo 
models, however, the solutions indeed form a coherent longitudinal structure over the 
whole latitudinal range, fulfilling the definition of an active longitude. Nevertheless, the 
paradigm of active longitudes has strongly influenced the discussion and analysis methods 
of solar activity tracers; very often bimodal distributions are the ones searched for, i.e. the 
distributions producing the strongest two-peaked distributions are normally ranked to 
be the most desirable. Even if not, the statistical significance of longitudinal clustering of 
any number of peaks in the histograms has been found to be very low. Furthermore, the 
best period calculated for the whole data set might not represent well both hemispheres, 
nor be the best if the data was divided into shorter intervals. This has led to the overall 
conclusion that the non-axisymmetric component, at least the one that would be rigidly 
rotating, of the solar dynamo must be very weak and apparently dynamically changing 
over time. 

It is well known that the solar convection zone is differentially rotating so that at 
the surface, the pole rotates slower than the equator. The next hypothesis adopted was 
to postulate, that the surface distribution of activity traces becomes affected by this 
rotation pattern. This essentially means that time series analysis methods allowing for 
one free parameter more to be fitted, namely the amount and direction of the surface 
differential rotation pattern, were developed. Again, using the maximal clustering as 
the measure of the statistical significance of the solution, it was claimed (Berdyugina & 
Usoskin, 2003; Usoskin et al., 2005; Berdyugina et al., 2006) that a persistent system of 
active longitudes, comprising roughly 10% of the total number of sunspots, was present 
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Figure 1. Photometric lightcurve of FK Comae Berenices for 1987-2004 analyzed with the Car- 
rier Fit method (Pelt et al., 2011; see also Hackman et al., 2012b). Number of model harmonics 
used is K = 3, number of modulation harmonics L — 20, modulation period of 9000 days, and 
carrier period P =2.4002466d. 

in the Sun, made visible by the 'passing' sunspot distribution on the surface experiencing 
differential rotation. The existence of this stroboscopic effect was put under doubt by the 
following investigations (Pelt et al. 2005; 2006), which showed that even the differentially 
rotating constructions produced longitudinal distributions that had a significance level 
comparable to a statistical fluctuation. On the other hand, it was shown (Pelt et al. 
2010) that the sunspot distribution, and therefore the magnetic field structure from 
which it originates, is indeed affected by differential rotation. Moreover, it was shown 
that nonparametric methods can be used to determine the coherence time of the non- 
axisymmetric structures; the average coherence time was estimated to be roughly 10-15 
Carrington rotations. 

As a conclusion, one might say that a weak non-axisymmetric component of the solar 
magnetic field exists, and is affected by the differential rotation of the solar convection 
zone. Due to its weakness and dynamical nature, it is very hard to quantify it even with 
the best statistical tools and despite of the long time span of the sunspot data. Most of 
the solar dynamo models work under the assumption of axisymmetry, which seems to 
be safe, taken what is said above. Here we note, however, that also non-axisymmetric 
solar dynamo models have been developed (e.g. Moss, 1999), where the excitation of a 
sub-dominant non-axisymmetric magnetic field modes has been seen. These models are 
kinematic, i.e. no back reaction of the magnetic field on the velocity field is taken into 
account. Therefore, these models do not properly address the dynamical significance of 
these modes. 

2.2. Rapidly rotating late-type stars 

Direct observational ways to detect and characterize spots and their longitudinal distri- 
butions on stellar surfaces include photometric lightcurve analysis, spectroscopic Doppler 
imaging, and spectropolarimetric Zeeman-Doppler imaging. Observational result that has 
had the strongest impact on the discussion of active longitudes came from the analysis 
of the photometry of FK Comae Berenices by Jetsu et al. (1993). By analyzing 25 years 
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of photometry, two active longitudes, roughly 180 degrees apart were found. The ac- 
tivity was seen to abruptly switch from one longitude to the other, the jump detected 
three times in the data used. This phenomenon was given the name 'flip-flop'; analysis of 
lightcurves of other rapidly rotating late-type stars (e.g. Berdyugina & Tuominen, 1998) 
soon established the fact that flip-flops are a common phenomenon related to the stellar 
cycle. Doppler imaging results, yielding the surface temperature distribution, confirmed 
the photometric results (e.g. Korhonen et al., 2004 and references therein), although with 
poorer time sampling. 

When these stars were followed up for a longer time span, it first became evident that 
the active longitude system is not stable in the rotational frame of reference of the stellar 
surface, i.e. the non-axisymmetric system either formed straight, inclined, railroads either 
upwards and downwards on the phase-time plots (e.g. Berdyugina & Tuominen, 1998), 
or even more complicated patterns of railroads with changing inclination (Jetsu et al., 
1999 for V 1794 Cyg), sinusoidal paths (Berdyugina et al. 2002 for LQ Hya; Berdyugina 
& Jarvinen, 2005 for AB Doradus) and so on. This revealed the other type of time- 
dependent behavior of the active longitudes - the non-axisymmetric system rotates with 
a different speed than the star itself, and this rotation period actually changes over time. 
This does not sound completely unfamiliar, as the solar hot spots are known to exhibit 
the same behavior. Secondly, the longer the stars were followed, the less periodic and 
regular the flip-flop events appeared. To illustrate the situation for FK Comae Berenices, 
we have re-analysed its photometry for the years 1987-2004 (Olah et al., 2006; Korhonen 
et al., 2007) with the Carrier Fit (CF) method (Pelt et al., 2011), see Fig. □ Here, we 
have used the well-known photometric rotation period of Pq ~ 2.4d as the carrier in our 
analysis, and computed the continuous longitudinal distribution of spots, appearing as 
dark regions in the plot, on the star. One of the 'well-behaved' flip-flop events reported by 
Jetsu et al. (1993) is visible during the first 3 years of the data. Since then, at least four 
other flip-flop type events have occurred, but the time intervals between them is clearly 
not regular. Moreover, railroad patterns (tracks tilted up- and downwards) appear in the 
data since 1995 indicating that the non-axisymmetric structure does no longer rotate 
with the photometric rotation period, but can lag behind or speed up with respect to 
it. This photometric analysis is in rough agreement with Doppler imaging results, see 
Korhonen et al. (2004). For a more detailed analysis of the phase jumps we refer to Olah 
et al. (2006) and Hackman et al. (2012b). 

The railroad pattern is especially clear and persistent in the case of the RS CVn binary 
system's primary giant component II Peg (see Fig. [2j Lindborg et al., 2013). Instead of 
two active longitudes, in this star we see the dominance of one dark region, with sporadic 
appearance of a secondary minimum roughly 180 degrees apart from the primary one (see 
Berdyugina et al., 1998, 1999; Lindborg et al., 2011; Hackman et al., 2012a for surface 
images). The dominance of the other active longitude could be thought to arise from 
the binarity, but the situation is far from being that simple. Namely, the CF analysis 
reveals a persistent downward trend in the phase-time diagram over most of the timespan, 
indicating that the spotted structure on a single active longitude is rotating faster than 
the star itself. This is a very peculiar finding, especially reflected upon the fact that in 
the system of close binary stars, the orbital period becomes synchronized. It is possible 
that the synchronization of the orbit has not yet fully occurred, but then the dominance 
of one active longitude due to binarity and synchronism appears invalid. Moreover, the 
systematic trend seems to break down towards the end of the data set. During this epoch, 
Doppler imaging (Hackman et al., 2012a) and Zeeman Doppler imaging (Kochukhov et 
al., 2012) suggest that the magnetic activity level of the star is attaining a minimum. 
Generally, the analysis of photometry is in good agreement with the Doppler imaging 



6 



M. J. Mantere, P. J. Kiipyla & J. Pelt 



Q- 




1990 1995 2000 2005 2010 

Time (years) 



Figure 2. Photometric lightcurve of II Peg for 1988-2010 analyzed with the Carrier Fit method 
(Pelt et al., 2011). Number of model harmonics used is K = 3, number of modulation harmonics 
L — 5, modulation period of 10000 days, and Carrier period Po~ 6.724333d. 

results (see e.g. Hackman et al., 2011, Lindborg et al., 2013), although the decreasing 
trend in the spot intensity is not visible in Fig. [5] due to the adopted normalization in 
our plot (each stripe is normalized with its cxtrema to enhance the drift pattern during 
the epoch of weaker spots) . 

Comparing photometry and Doppler images to magnetic maps (Kochukhov et al., 
2012) brings further complications into the interpretation, although it has to be kept in 
mind that magnetic maps exist only for the lower state activity epoch - magnetic maps 
clearly show two active longitudes of different polarities virtually at all times the star 
has been imaged, which undergo polarity changes without any signs of flip-flop occurring 
with the other methods. Only exception to this are the most recent observing seasons 
(2008-2010), during which a polarity change might be associated with the abrupt-looking 
phase behavior seen also in Fig. [5J 

3. Explanations in the framework of turbulent dynamo theory 

To understand the gradually building up observational picture, theory should address 
at least the following two fundamental issues: Firstly, the models should be capable 
of reproducing the transition from oscillatory, mostly axisymmetric dynamo in the afl 
regime (differential rotation dominating over the effect of convective turbulence) into the 
dominant ly non- axisymmetric dynamo in the pure a 2 regime (field generation presumably 
taking place with negligible differential rotation) when the rotation rate of the star is 
increased. The usage of 'presumably' here refers to the fact that we do not actually 
yet observationally know the internal rotation profiles for other stars than the Sun; the 
only information comes from modeling (see e.g. Kitchatinov & Riidiger, 1999). Secondly, 
what is the relation of the large-scale dynamo-generated magnetic field and the sun- and 
starspots themselves, and how does one form spots from the underlying dynamo field? 

Both of these issues have been addressed from the point of view of a turbulent dynamo, 
which viewpoint is adopted in this paper; the competing theory of the transport of rising 
flux tubes is discussed in length by Choudhuri (2013) in these proceedings. According 
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Figure 3. Radial velocity U r from three simulations of rapidly rotating turbulent convection 
with density contrasts of 5 (left), 30 (middle), and 100 (right) in spherical wedge geometry. 
Adapted from iKapyla et all l|2011br i. 

to the theory of a turbulent dynamo (for an extensive review, see Brandenburg & Sub- 
ramanian, 2005), the dynamo instability itself is affected by the presence of turbulence 
so that a collective inductive action plus an enhanced diffusion is caused by turbulence 
effects; as the convection zone is turbulent throughout, these effects arc distributed over 
the total extent of it - sometimes the term distributed dynamo is used to denote the 
turbulent dynamo process. Turbulence effects can also lead to certain types of instabil- 
ities, that have been shown to result in the formation of magnetic field concentrations 
or global-scale temperature anomalies. On the other hand, structure formation in con- 
vectively turbulent flow can occur via other, less specified, routes, which provides us a 
starting point for the theoretical discussion. 

3.1. Spontaneous formation of magnetic field concentrations in turbulent convection 

simulations 

Direct numerical solutions of magnetized convection have reached a level of sophistication 
where the formation of flux concentrations can be self-consistently modeled. Global and 
semi-global simulations of rapidly rotating convection exhibit non-ax isymmetric convec- 



tion at low latitudes when the magnetic fie ld is weak or absent (e.g. iBrown et al.l 12008 



Kapyla, et al.ll2011bt lAugustson et al.ll2012l see the leftmost panel of Fig. [3]). This be- 



havior is most prominent near the onset of convection. Furthermore, in creasing den- 



sity st ratification also seems to suppress the non-axisymmetric modes (jKapvla et al 



2011bl compare the middle and rightmost panel of Fig. [3] to the leftmost one). How- 



ever, as the magnetic fields due to dynamo action become dynamically important, the 
non-axisymmetric convection pattern disappears and the large-scale magnetic fields are 



also mostly axisymmetric (e.g. lBrown et al.ll201ll) : see, however, Miesch et al. (2013) in 
these proceedings. 

Another line of study concerns local simulations where some sort of magnetic field is 
imposed or advected into the system. In these cases turbulent convection can rearrange 
the magnetic field self-consistently into concentrations in contrast to the sunspot simu- 
lations of Rcmpcl et al. (2009) where the spots are the result of an imposed magnetic 
field structure. The outcomes of these models are either large-scale magnetic structures 
(e.g. Tao et al., 1998), or can resemble pores (Kitiashvili et al., 2010) or even bipolar 
regions in the Sun (Stein & Nordlund, 2012). The exact mechanism responsible for the 
flux concentrations in these cases is not yet certain. 

3.2. NEMPI 

A possible mechanism for generating flux concentrations that ultimately lead to the 
formation of sunspot and starspot s is due to a negative contribution of turbulence t o 
the effective magnetic pressure (e.g. lKleeorin et al. Hl99rtlRogachevskii fc Kleeorirj|2007t) . 
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Figure 4. Large-scale vortices in rapidly rotating turbulent convection in Cartesian geometry. 
Left and middle panels show the temperature and vertical velocity, respectively, near the surface 
of the convectively unstable region. The right panel shows a cut of the temperature thro ugh the 
warm spot in the lower left corner of the left panel. Adapted from iMantere et alJ (|201lh . 

a.k.a. the negative effective magnetic pressure instability (NEMPI). This effect can form 
flux concentrations even from uniform, sub-equipartition, magnetic fields. 

Mean-field models and direct numeric al simulations of forced turbule nce have estab- 
lished that the negative effect exists (e.g. Brandenburg et aLll2010l 12012 ). and the insta- 



bility itself has been detec ted in direct numerical simulations (DNS) of forced turbulence 



( Brandenburg et aill2011 ) . A neg ative effect on magnetic pressure has also been found for 



convection (jKapyla et al.ll2012bf ) but no NEMPI has been detected so far (see, however, 
Kapyla et al. (2013) in these proceedings). 

3.3. Large-scale hydrodynamic vortex instability 

The apparent uncorrelation of magnetic fields and temperature anomalies in Zeeman- 
Doppler maps of stars is a puzzling feature, which may be attributable to the limited 
spatial resolution of the current ZDI maps (see discussion in Kochukhov et al., 2012), 
or explained with a new mechanism by which starspots can be generated without a 
correlation of magnetic fields and temperature. The latter possibility has been recently 
considered on the basi s of results from loca l hydrodynamic simul ations of rapidly rotating 
turbulent convection (Ka pyla, et al. 2011a ; Mantere et al.|[20Tlh . It turns out that when 



certain threshold values of the the Reynolds and Coriolis numbers, describing the effects 
of molecular viscosity and rotation on the flow, are exceeded, large-scale vortices appear 
in the system. When rotation is gradually increased, cool anti-clockwise vortices appear 
first. For more rapid rotation also warm clockwise cyclones appear (see Fig. [4j. The 
temperature anomaly in the cyclones is at least ten per cent of the surface temperature 
which similar to the spots seen in Doppler images of rapidly rotating stars. Similar 
vortices have not yet been found in global or semi-global simulations. A probable reason 
is too low resolution achievable at the moment. Furthermore, t he vortices also seem 



to disappear in simulations where a magnetic dynamo is present (jKapyla et al.ll2012al ). 
Currently it is unclear whether the cyclones can coexist with strong magnetic fields. 

3.4. Non-axisymmetric mean-field dynamos 

Modeling of the solar dynamo is at the verge of becoming accessible by DNS (Ghizaru et 
al. 2010), the first solar- type solutions having been obtained from such models (Kapyla 
et al., 2012c). Such computations, however, are computationally expensive, and mean- 
field modeling still has clear advantages in this repspect. The solar mean-field models 
usually work under the assumption of axisymmetry, reducing the problem into two spatial 
dimensions, whereas it is obvious from observations that this assumption is not valid for 
the more rapidly rotating stars with dominantly non-axisymmetric surface temperature 
and magnetic field configurations. Non-axisymmetric models necessarily solve for the 
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azimuthal dependence of the physical quantities, either on a finite-difference or -volume 
grid, or using spherical harmonic decomposition and Fourier transform. 

The non-axisymmetric solutions to the kinematic a 2 equations in simple setups were 
computed already decades ago (see e.g. Radler, 1975; Krause & Radler, 1980), providing 
more or less all the essential properties of the solutions: when rotation rate is increased 
high enough, the preferentially excited modes are the non-axisymmetric ones, typically 
representing quadrupolar symmetry (SI). The solutions are usually non-oscillatory, but 
still time-dependent, as the modes migrate either east- or westward. The nonlinear, some- 
times even dynamical, models developed since the 1990's (see e.g. Moss et al. 1995; Kiiker 
& Riidiger 1999) have shown that also antisymmetric dipolar solutions (Al) can appear 
as the nonlinearly stable solutions (Tuominen et al. 2002). In this type of a solution, 
two high-latitude spots of opposite polarities in the same hemisphere are generated, 
in agreement with what is observed e.g. in II Peg (e.g. Kochukhov et al., 2012). De- 
spite of the azimuthal dynamo wave generated by the east- or westward migration of 
the non-axisymmetric modes, the solutions reported have been non-oscillatory, therefore 
providing no explanation for the flip-flop type phase jumps. 

In search for oscillatory solutions, one approach has been to include small amounts 
of differential rotation into the models (see e.g. Elstner & Korhonen, 2005), essentially 
meaning that the dynamo models are of a 2 Q type, where differential rotation is present, 
but is subdominant in the magnetig field generation. In this way, mixed solutions with 
an oscillatory axisymmetric mode and a stable non-axisymmetric one provide flip-flop 
type oscillatory behavior. On the other hand, it has been known since a long time that 
pure a 2 dynamos do have oscillatory solutions (see e.g. Baryshnikova & Shukurov, 1987; 
Brauer & Radler, 1987). In fact, only in the case if the a effect is purely homogeneous and 
the system unbounded, the solutions are purely non-oscillatory. Recently, oscillatory a 2 
solutions have been found in the spherical forced turbulence simulations of Mitra et al. 
(2010) and in local Cartesian simulations of turbulent convection (Kapyla et al., 2012a). 
In the next subsection, we illustrate the possibility of oscillatory non-axisymmetric a 2 
solutions in a simple mean-field dynamo model. 



3.4.1. Simple a 2 dynamo model 

We solve the mean-field induction equation in three dimensions using spherical polar 
coordinates (r, 9, 4>). The domain is a wedge spanning the latitude range (6°, 174°), radius 
0.6, 1.0 and the full azimuthal extent. The equation is solved with the Pencil CODE0 
using a mesh 32 x 64 x 128; the spatial discretization in spherical polar coordinates does 
not allow us to explore the full latitudinal extent, as the time step decreases strongly close 
to the pole. The a effect has a uniform radial profile, and cos£? dependence, obtaining 
maxima at the poles. The boundary conditions applied for the magnetic field, in terms 
of the vector potential that is solved for, are those of a perfect conductor at the lower 
radial and latitudinal boundaries, and normal-field conditions are applied at the top 
radial boundary: 

(r = n), (3.1) 
(r = r 3 ), (3.2) 
(0 = 0i,02). (3.3) 



dAr 
dr 

A r 



Ad = A, 



0,- 



dA g 





Ag dA d 



dr 



OAg 

88 



A 6 = 



t http:/ /pencil-code. googlecode.com/ 
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Figure 5. Radial magnetic field at the surface for six different times, 
t/T=(4.17,4.28,4.39,4.50,4.61,4.72) of the oscillation cycle of a solution with C a =5 x 10 3 . 

Azimuthal boundaries are treated as periodic. As a nonlinearity, a simple algebraic a- 
quenching is used. 

We perform calculations with varying strength of the a-effect, spatially constant tur- 
bulent diffusion, and time being measured in diffusion times, r = d 2 /r], where d is the 
depth of the convection zone. The strength of the a effect is measured by the quantity 
C a = aoR/r], where ao is varied. The critical C a above which growing solutions are 
obtained is roughly 50. A regime of axisymmetric oscillatory solutions of antisymmetric 
type is seen up to C a « 10 3 . With larger values, the excited solutions become dom- 
inantly non-axisymmetric, and show eastward migration. The solutions are otherwise 
non-oscillatory. When a effect is ramped up further, the solutions change into oscilla- 
tory ones. Time evolution of the radial field at the surface during one cycle of such a 
calculation is shown in Fig. [5J and magnetic energy densities for different magnetic field 
components in Fig. [5] The oscillation seen in the magnetic field energy densities is clearly 
related to the interchange of the magnetic field strength of spots with different polarities. 
Interestingly, the migration period of eastward drift equals the cycle length. 



4. Conclusions 

In this paper we have discussed observations and theory of longitudinal concentrations 
of solar and stellar activity tracers. While the solar magnetic field is dominantly axisym- 
metric, non-axisymmetric modes dominate the magnetic fields of rapid rotators. The 
greatest challenges in the observational frontier include understanding the uncorrelation 
of temperature anomalies from magnetic field strength in ZDI, and determining whether 
and how the flip-flops and phase jumps relate to stellar activity cycles. Theoretically 
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Figure 6. Left: snapshot of the radial magnetic field at the surface for i/r=4.39 in spherical 
projection. Right: magnetic energy density of the different components (total -thick solid, radial 
- thin solid, latitudinal - dotted, azimuthal - dashed). 

the most challenging task consists of explaining how the global dynamo-generated field 
transforms itself into sun- and starspots, and whether the mechanism is the same in all 
stars independent of the rotation rate. 
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